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Abstract This work investigates the muon capture reactions H(/i , f^Jnn and 
3 He(/x _ ,^) 3 H and the contribution to their total capture rates arising from the ax- 
ial two-body currents obtained imposing the partially-conserved-axial-current (PCAC) 
X^/y ' hypothesis. The initial and final A = 2 and 3 nuclear wave functions are obtained from 

the Argonne two-nucleon potential, in combination with the Urbana IX three- 
nucleon potential in the case of A = 3. The weak current consists of vector and axial 
components derived in chiral effective field theory. The low-energy constant entering 
the vector (axial) component is determined by reproducting the isovector combination 
of the trinucleon magnetic moment (Gamow- Teller matrix element of tritium beta- 
decay). The total capture rates are 393.1(8) s _1 for A = 2 and 1488(9) s _1 for A = 3, 
■ where the uncertainties arise from the adopted fitting procedure. 
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1 Introduction 



o 

■ There is a significant body of experimental and theoretical work on muon captures 

in light nuclei, motivated by the fact that the theoretical framework used to study 
these reactions is the same as that used for weak capture reactions of astrophysical 
interest, not accessible experimentally. Muon captures, whose rates can be measured, 
can therefore provide a valuable test of this theoretical framework pQ. 
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Very recently [2], the muon capture reactions 2 H(/i~, ^ t )nn and 3 He(/i~, ^) 3 H 
have been studied simultaneously in a consistent framework. In particular, the initial 
and final A = 2 and 3 nuclear wave functions have been obtained from the Argonne 
«18 (AV18) [3j or the chiral N 3 LO (N3LO) [3] two-nucleon potential, in combination 
with, respectively, the Urbana IX (UIX) [5] or chiral N 2 LO (N2LO) [6] three-nucleon 
potentials in the case of A = 3. The weak current consists of polar- and axial- vector 
components. The former are related to the isovector piece of the electromagnetic cur- 
rent via the conserved-vector-current (CVC) hypothesis. These and the axial current 
have been derived within two different frameworks, the standard nuclear physics ap- 
proach (SNPA), and chiral effective field theory. The first one goes beyond the impulse 
approximation, by including meson-exchange current contributions and terms arising 
from the excitation of Z\-isobar degrees of freedom. The second approach includes 
two-body contributions derived in heavy-baryon chiral perturbation theory within a 
systematic expansion, up to N 3 LO [7]|8]. To be noticed that, since the transition oper- 
ator matrix elements are calculated using phenomenological wave functions, it should 
be viewed as a hybrid chiral effective field theory approach (EFT*). Both SNPA and 
EFT* frameworks have been used in studies of weak pp and hep capture reactions in 
the energy regime relevant to astrophysics [71l9lll0]. The only parameter in the SNPA 
nuclear weak current model is present in the axial current (the iV-to-Zi axial coupling 
constant) and is determined by fitting the experimental Gamow- Teller matrix element 
in tritium /3-decay (GT EXP ). The SNPA weak vector current, related to the isovec- 
tor electromagnetic current via CVC, reproduces the trinucleon magnetic moments to 
better than 1 % [2]. In the case of EFT*, three low-energy constants (LECs) appear: 
one in the axial-vector component, and two in the electromagnetic current. Of these, 
only one is relevant to the weak vector current, since the other appears in front of 
an isoscalar operator. The corresponding coupling constants are parameters fixed to 
reproduce, respectively, GT EXP and A — 3 magnetic moments. To be noticed that 
the EFT* currents are obtained performing the Fourier transform from momentum- to 
coordinate-space with a Gaussian regulator characterized by a cutoff A, varied between 
500 and 800 MeV. The total capture rates have been found to be 392.0(2.3) s^ 1 for 
A = 2 and 1484(13) s _1 for A = 3. The spread accounts for the model dependence, i.e., 
the dependence on the input Hamiltonian model, the model for the nuclear transition 
operator, and, in the EFT* calculation, the cutoff sensitivity. This weak model depen- 
dence is a consequence of the procedure adopted to constrain the weak current. These 
results are in very good agreement with the experimental data, in particular with the 
very accurate measurement of Ref. [11] for the total rate in muon capture on 3 He. 

The muon capture on deuteron has been studied also using the SNPA and the EFT* 
framework in Ref. [12] . The SNPA retains two-body meson- exchange currents derived 
from the hard pion chiral Lagrangians of the N Air puia\ system and are significantly 
different from those of Ref. [2]- On the other hand, the EFT* currents are similar to 
those of Ref. [5], but two differences need to be remarked: (i) the LEC appearing in 
the axial- vector component (dj{) is not fixed to reproduce GT EXP , rather the doublet 
capture rate calculated in SNPA; (ii) a term is added to the leading axial two-body 
currents, in order to satisfy the partially-conserved-axial-current (PCAC) hypothesis, 
as constructed in Ref. Q3] (called there, and from now on, potential current). The 
calculated SNPA values for the total capture rate are in the range of 416-430 s _1 , de- 
pending on the potential model used, resulting in a model dependence much larger than 
in Ref. 2 . It is also argued that "omitting the potential current causes an enhancement 
of the doublet transition rate ^1/2 by ~ 1%" [12] , 
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In the present work we repeat the calculation of Ref. [2], in the EFT* approach, 
adding the potential currents as in Ref. [12| . We restrict our calculation to the AV18 
and AV18/UIX potential models, and to a cutoff value of 600 MeV, since, as shown 
in Ref. [2], the dependence on these inputs is less than 1 %. We fit the d,R coefficient 
to GT EXP , and consistently calculate the total rates for muon capture on deuteron 
and 3 He. The comparison with the results of Ref. 2 will give and indication of how 
significant are the potential current contributions for these muon captures. 

The paper is organized as follows: in Sec. [2] the theoretical formalism used in the 
calculation is briefly reviewed. In Sec. [3] the EFT* model for the weak current is de- 
scribed, with the addition of the potential currents. In Sec.0] the results are presented 
and discussed, and some concluding remarks are given. 



2 Theoretical formalism 

We briefly review the formalism used in the calculation for the muon capture processes, 
discussed at length in Refs. [2, 14 . The muon capture on deuteron and 3 He is induced 
by the weak interaction Hamiltonian IS], Hyy = ^K- J dxZ cr (x)j' T (x), where Gy is 

the Fermi coupling constant, G v = l. 14939 xl0~ 5 GeV~ 2 16 , and l a and f are the 
leptonic and hadronic current densities, respectively. The transition amplitude can be 
written as 

Tw(f,fz;s 1 ,s 2 ,h l/ ) = {nn,s 1 ,s 2 ;v, h v \H W \ (p, d); /, f z ) 
~ %CX1^M.^I/M^(^ ^)<^i S 2(nn)|r(q)l^(s d )) , (1) 

for muon capture on deuteron, p being the nn relative momentum, and [14] 
T W (f, h;s' 3 , h v ) = < 3 H, 4 v, h u \H w \ (/i, 3 He); /, f z ) 

- E4^'^ S3|//2);ct( ^' ^)(*'3 H (4)b" T (q)I^Ho(s3)> , (2) 

for muon capture on 3 He. In order to account for the hyperfine structure in the initial 
system, the muon and deuteron or 3 He spins are coupled to states with total spin 
/ = 1/2 or 3/2 in the deuteron case, and f — or 1 in the 3 He case. In Eqs. |T| 
and ([2]) we have defined with {h y ) the muon spin (muon neutrino helicity). The 
Fourier transform of the nuclear weak current has been introduced as 

j CT (q) = y'dxe i< l x r(x) = (p(q),j(q)) , (3) 

with the leptonic momentum transfer q defined as q — ku — — — k&/, ku and kjy 
being the muon and muon neutrino momenta. The function tpfl has been introduced 
to take into account the initial bound state of the muon in the atom and the charge 
distribution of the nucleus. It is typically approximated as [15] l^ill 2 — ( a ft^) f or 

muon capture on deuteron, and Q3] \4>il \ 2 = ^ 3Hc ^ for muon capture on 3 He, 

where a is the fine structure constant (a = 1/137), fj,^ and /i^3 He are the reduced 
masses of the (/x, d) and (/x, 3 He) systems, and the factor 1Z approximately accounts for 
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the finite extent of the nuclear charge distribution |15| and is taken to be 0.98, as in 
Ref. 0. 

In the case of muon capture on deuteron, the final state wave function $p,sis 2 (nn) is 
expanded in partial waves, and the calculation is restricted to total angular momentum 
J < 2 and orbital angular momentum L < 3, i.e., in a spectroscopic notation, to 
^o, 3 P , 3 Pl, 3 P2~ 3 F 2 and 1 D 2 . Standard techniques [T0lfT5] are now used to carry 
out the multipole expansion of the weak charge, p(q), and current, j(q), operators. 
Details of the calculation can be found in Ref. [2]. Here we only note that all the 
contributing multipole operators selected by parity and angular momentum selection 
rules are included, as explained in Ref. [2]. 

The total capture rate for the two reactions under consideration is then defined as 

dr = 2ttS(AE)\T w \ 2 x (phase space) , (4) 

where 8{AE) is the energy-conserving 5-function, and the phase space is dp dk„/(27r) 6 
for muon capture on deuteron and just dk„/(27r) 3 for muon capture on 3 He. The 
following notation has been introduced: (i) for muon capture on deuteron 

^^=2/^1 ^2 Dl T w(/,/*;si, S 2,MI 2 . (5) 

and the initial hyperfine state has been fixed to be / = 1/2; (ii) for muon capture on 
3 He 

\ T W? = \ EE^^'^^MI 2 , (6) 

s' 3 h v ffz 

and the factor 1/4 follows from assigning the same probability to all different hyperfine 
states. 

After carrying out the spin sums, the differential rate for muon capture on deuteron 
(dF D /dp) and the total rate for muon capture on 3 He (To) are easily obtained, and 
their expressions can be found in Ref. [2]. In order to obtain the total rate F D for 
muon capture on deuteron, dr D /dp is plotted versus p and numerically integrated. 

Bound and continuum wave functions for both two- and three-nucleon systems 
entering in Eqs. JT]) and @ are obtained with the hyperspherical-harmonics (HH) 
expansion method. This method, as implemented in the case of A = 3 systems, has 
been reviewed in considerable detail in a series of recent publications 17,18:19 . We 
have used the same method in the context of A = 2 systems, for which of course wave 
functions could have been obtained by direct solution of the Schrodinger equation. A 
detailed discussion for the A — 2 wave functions is given in Ref. [2]. 

3 The nuclear weak current operator 

The chiral effective field theory weak current transition operator is taken from Refs. [7] 
and [8], as reviewed in Ref. [2]. It is derived in covariant perturbation theory based 
on the heavy-baryon formulation of chiral Lagrangians by retaining corrections up to 
N 3 LO. The one-body operators are those listed in Eqs. (17) of Ref. [7] and (4.13)-(4.14) 
of Ref. [2]. The vector charge and axial current operators retain terms up to 1/m 2 , while 
the axial charge and vector current operators retain terms up to 1/m 3 , m being the 
nucleon mass. Both 1/m 2 and 1/m 3 contributions arise when the non-relativistic re- 
duction of the single-nucleon covariant current is pushed to next-to-leading order. The 
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two-body vector currents are obtained from the two-body electromagnetic currents via 
CVC. These are decomposed into four terms [8J: the soft one-pion-exchange (f7r) term, 
vertex corrections to the one-pion exchange (ItvC), the two-pion exchange (2n), and 
a contact-term contribution. Their explicit expressions can be found in Ref. [HJ. All 
the 1-7T, InC and 2tv contributions contain low-energy constants estimated using reso- 
nance saturation arguments, and Yukawa functions obtained by performing the Fourier 
transform from momentum- to coordinate-space with a Gaussian regulator character- 
ized by a cutoff A. Here, as discussed above, we fixed the value of A = 600 MeV. The 
contact-term electromagnetic contribution is given as sum of two terms, isoscalar and 
isovector, each one with a coefficient in front (g^s and g^v) fixed to reproduce the 
experimental values of triton and 3 He magnetic moments. For the AV18/UIX Hamil- 
tonian model, with A = 600 MeV, (745 = 0.55(f) and g^y = 0.793(6), the error being 
due to numerics [2]- Note that only the isovector contribution is of interest here, but 
anyway it turns out to be negligible. 

The two-body axial current operator consists of two contributions: a one-pion ex- 
change term and a (non-derivative) two-nucleon contact-term. The explicit expression 
for the contact term can be found in Ref. [JJ. Here we review the one-pion exchange 
term, since we add, in accordance with Ref. |f 3| . the potential current contributions. 
Therefore, in momentum-space, the one-pion exchange term reads: 



I ( 1 - ^i) + ( i + -4 ) ^ x k. 



1 + c 6 + g\ 

+( -. — x q 



- 2£t/ (q + io-i x fa + p^)) 1 7 ' k * + i j . (7) 



J 



where kjj = p[ ■ — Pi.j, with pj j and p^ ■ being the initial and final single nucleon 
momenta, q = + k, , g^ = f .2654 is the axial-vector coupling constant, and f n = 93 
MeV is the pion decay constant. The values used for the coupling constants £3, £4, 
and cq, as obtained from nN data, are £3 = —3.66, £4 = 2.11 and cq = 5.83 [7]. The 
terms proportional to g\ in Eq. are the potential currents. They are the same 
as in Eq. (21) of Ref. [13] or Eq. (A.17) of Ref. [12]. The low-energy constant d R , 
determining the strength of the contact-term two-body axial contribution, has been 
fixed by reproducing GT EXP , finding dji — 1.54(8). This value should be compared 
with the corresponding one given in Ref. [2] (see Table V), <1r — 1.75(8). The difference 
between these two values of ~ 13 % is due to the presence of the potential currents 
and is comparable with that of Ref. [12] . 



4 Results 

We present in Table [T] the results for the total rates of muon capture on deuteron, 
in the doublet hyperfine state (F D ), and on 3 He (Jo). The deuteron, nn, 3 He and 
3 H wave functions have been calculated with the AV18 [3] two- and, when necessary, 
UIX [5] three-nucleon interactions. The model for the nuclear weak transition operator 
has been presented in Sec. [3] We compare our results with those of Ref. [2], obtained 
with the same Hamiltonian model and cutoff A, but without the two-body potential 
currents elaborated in Ref. [13] and discussed in Sec. [3] Note that here dfj = 1.54(8), 
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Table 1 Total rate for muon capture on deuteron and 3 He, in s _1 . In the A = 2 case, the 
different partial wave contributions are indicated. The numbers among parentheses indicate 
the theoretical uncertainties arising from the adopted fitting procedures. Such uncertainty is 
not indicated when less than 0.1 s _1 . The AV18 and AV18/UIX interactions have been used 



to calculate the A = 2 and A = 3 wave 
also listed. 


functions 


The 


corrct 


ponding results of Ref. 




1 S 3 Po 3 Pi 


3 P 2 


X D 2 


3 P 2 




To 


Present work 
Ref. [2] 


250.1(8) 20.2 46.1 
250.0(8) 19.8 46.3 


71.3 
71.1 


4.5 
4.5 


0.9 
0.9 


393.1(8) 
392.6(8) 


1488(9) 
1488(9) 



while in Ref. [5] d,R — 1.75(8). From inspection of the table, we conclude that the two 
calculations are in remarkable agreement with each other: the differences in Po an< l 
r D arc < 0.1 %, well below the theoretical uncertainties. The largest difference, of the 
order of 2 %, is in the 3 Prj partial wave contribution to r D . However, when all the 
partial wave contributions are summed up, the difference in F D returns well below the 
1 % level. 

In conclusion, we have studied the potential currents dictated by PCAC, as elab- 
orated in Refs. [121113] . and we have found that their contributions to the total rates 
of muon capture on deuteron and 3 He are tiny. This result is a consequence of the 
procedure adopted to constrain the weak current. Finally, we expect that the potential 
currents will give tiny contributions also in weak capture reactions of astrophysical 
interest and in those processes whose momentum transfer is small. 
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